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ABSTRACT 

This paper uses the First XMM-Newton Serendipitous Source Catalog compiled by 
£ — ' the XMM-Newton Science Center to identify \ow-z X-ray selected normal galaxy 

| candidates. Our sample covers a total area of « 6deg 2 to the 0.5-2 keV limit 

£NJ . w 10 -15 ergs -1 cm -2 . A total of 23 sources are selected on the basis of low X-ray-to- 

optical flux ratio log fx/ f opt < — 2, soft X-ray spectral properties and optical spectra, 
C> when available, consistent with stellar than AGN processes. This sample is combined 

with similarly selected systems from the Needles in the Haystack Survey (Georgan- 
topoulos et al. 2005) to provide a total of 46 unique z 55 0.2 X-ray detected normal 
r"| . galaxies, the largest low-z sample yet available. This is first used to constrain the nor- 

rnal galaxy logiV — logS at bright fluxes (10 -15 — 10 -13 erg s -1 cm -2 ). We estimate a 
slope of —1.46 ± 0.13 for the cumulative number counts consistent with the euclidean 
prediction. We further combine our sample with 23 local (z ^ 0.2) galaxies from the 
Chandra Deep Field North and South surveys to construct the local X-ray luminosity 
, function of normal galaxies. A Schcchter form provides a good fit to the data with a 

break at logL* = 41.02^° ergs -1 and a slope of a — —1.76 ± 0.10. Finally, for the 
sample of 46 systems we explore the association between X-ray luminosity and host 
galaxy properties, such as star-formation rate and stellar mass. We find that the Lx 
5^ ■ of the emission-line systems correlates with Ha luminosity and 1.4 GHz radio power, 

both providing an estimate of the current star-formation rate. In the case of early type 
galaxies with absorption line optical spectra we use the if-band as proxy to stellar 
mass and find a correlation of the form Lx oc L]^ 5 . This is flatter than the Lx — Lb 
relation for local ellipticals. This may be due to either Lk providing a better proxy to 
galaxy mass or selection effects biasing our sample against very luminous early-type 
galaxies, Lx > 10 42 ergs -1 . 
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1 INTRODUCTION 

The launch of the XMM-Newton and the Chandra missions 
has opened new opportunities in the study of the X-ray 
properties of normal galaxies. Apart from targeted observa- 
tions of local systems, surveys performed by these telescopes 
have allowed, for the first time, study of normal galaxies at 
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X-ray wavelengths outside the nearby Universe, thus open- 
ing the way for evolutionary investigations. Stacking analy- 
sis studies using these missions for example, have provided 
constraints on the mean X-ray properties {Lx , hardness ra- 
tios, Lx I Lb) of star- forming systems over the redshift range 
z = 0.1 - 3 (Hornschemeier et al. 2002; Nandra et al. 2002; 
Georgakakis et al. 2003a, b; Laird et al. 2004). The stacking 
analysis results from the above independent studies when 
combined together are consistent with luminosity evolution 
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of the form « (1 + z) 3 at least tozfs 1.5 (e.g. Georgakakis 
ct al. 2003a, b). 

The X-ray stacking analysis provides information on the 
mean properties of systems that are too X-ray faint to be in- 
dividually detected. A major breakthrough however, of the 
XMM-Newton and the Chandra missions has been the de- 
tection of normal galaxies at low and moderate redshifts. 
Bauer et al. (2002) and Alexander et al. (2002) studied the 
radio and mid-infrared properties of X-ray sources in the 
Chandra Deep Field (CDF) North and argue that systems 
with X-ray emission dominated by starburst activity are de- 
tected in this field out to z « 1. Hornschemeier et al. (2003) 
demonstrated that more quiescent normal galaxies with X- 
ray-to-optical flux ratio log fx /fopt < —2 and a median 
redshift z w 0.3 not only constitute a non-negligible com- 
ponent of the X-ray source population in the CDF-North 
but are also likely to outnumber AGNs at fluxes below 
/x(0.5-2keV) « 10~ 17 ergs" 1 cm" 2 . Norman et al. (2004) 
extended the Hornschemeier et al. (2003) study and identi- 
fied over 100 starburst/quiescent galaxy candidates to z w 1 
in the combined CDF-North and South albeit with optical 
spectroscopy limited to a fraction of them. These authors es- 
timate the X-ray luminosity function of normal galaxy can- 
didates and find evidence for X-ray evolution of the form 
tx (1 + z) 3 to 2 ?s 1, i.e. similar to that inferred from other 
wavelengths (e.g. Hopkins 2004 and references therein). 

Parallel to the studies above focusing on distant sources 
(z <; 0.3), there have also been efforts to compile X-ray se- 
lected normal galaxy samples at low redshifts z < 0.2. This 
is essential to complement the deep CDF studies and to 
provide a comparison sample in the nearby Universe. Geor- 
gakakis et al. (2004) and Georgantopoulos, Georgakakis & 
Koulouridis (2005) used public XMM-Newton data overlap- 
ping with the Sloan Digital Sky Survey (SDSS; Stoughton 
et al. 2002; York et al. 2000) to identify normal galaxies 
at a median redshift of about 0.07 over an lldeg 2 area. 
This dataset, the Needles in the Haystack Survey, has been 
used to constrain the galaxy log N — log S at bright fluxes, 
estimate the X-ray luminosity function of absorption and 
emission line galaxies separately and to assess their con- 
tribution to the diffuse X-ray background under different 
evolution scenarios. A similar study using the HRI detec- 
tor onboard ROSAT has been performed by Tajer et al. 
(2005) aiming to determine the number counts of low X- 
ray-to-optical flux ratio systems (log fx/ fopt < — 1) at very 
bright fluxes (10~ 14 — 10~ 12 erg s _1 cm -2 ). More recently, 
Hornschemeier et al. (2005) explored the X-ray properties 
of 2 ps 0.1 galaxies identified in public Chandra pointings 
overlapping with the SDSS. They use this sample to inves- 
tigate the link between X-ray luminosity, host galaxy star- 
formation rate (SFR) and stellar mass. These relations are 
consistent with a picture where the X-ray emission arises in 
a combination of low- and high-mass X-ray binaries as well 
as hot gas interstellar medium. 

Despite significance progress, the number of X-ray se- 
lected normal galaxies at low-2 (e.g. z ^ 0.2) remains 
small hampering statistical studies. In this paper we fur- 
ther expand existing samples using as a resource the 1st 
XMM-Newton Serendipitous Source Catalogue produced by 
the XMM-Newton Science Center. Our aim is to demon- 
strate the power of this unique and continuously expanding 
database and how it can be exploited to advance our under- 



standing of the properties of X-ray selected normal galax- 
ies. Among others, such a sample can provide a tight anchor 
point for comparison with deeper surveys probing on aver- 
age higher-2, calibrate the relation between SFR diagnos- 
tics and X-ray luminosity and improve our understanding of 
the association between host galaxy properties (e.g. stellar 
mass) and X-ray emission. Throughout this paper we adopt 
H = 70kms~ 1 Mpc~\ Q. M = 0.3 and ft A = 0.7. 

2 SAMPLE SELECTION AND OBSERVATIONS 
2.1 The X-ray data 

In this paper we use the First XMM-Newton Serendipitous 
Source Catalogue (1 XMM; version 1.0.1) constructed by the 
XMM-Newton Survey Science Centre (SSCT) and released 
on 2003 April 7th. This is a huge database comprising source 
detections drawn from 585 XMM-Newton EPIC observations 
made between 2000 March 1 and 2002 May 5. 

After the initial reduction of the raw data and the fil- 
tering of the high particle background periods the resulting 
event files are used to construct images in 5 energy bands: 
0.2-0.5, 0.5-2.0, 2.0-4.5, 4.5-7.5 and 7.5-12.0 keV. The ini- 
tial source detection is performed simultaneously on the 5 
energy bands above using the eboxdetect task of SAS. The 
detected sources are passed on to the SAS task emldetect 
to assess their reliability and to determine various source 
parameters by fitting the instrumental point spread func- 
tion (PSF). The source position was fit simultaneously on 
all input images, whereas the source count rates are left to 
vary independently in each image. Emldetect estimates 
count rates corrected for vignetting, losses due to inter-chip 
gaps, bad pixels/columns, events arriving during readout 
times and the extended PSF. Parallel to the source count 
rate the emldetect task also outputs the detection like- 
lihood for each source providing an estimate of the proba- 
bility of spurious detection. We convert count rates to flux 
assuming a power-law spectral energy distribution with in- 
dex r = 1.8 and Galactic absorption appropriate for each 
individual XMM-Newton field. The index V — 1.8 is consis- 
tent with the median hardness ratio, HR « —0.6, of late- 
type galaxies in our final sample (see section 3). Using a 
Raymond-Smith model with a temperature of 1 keV instead 
of a power-law, would overestimate the flux in the 0.5-2 keV 
band by only about 4 per cent. 

In this paper we concentrate on XMM-Newton fields 
that have the EPIC (European Photon Imaging Camera; 
Striidcr et al. 2001; Turner et al. 2001) cameras as the 
prime instrument with the PN detector operated in full- 
frame mode and exposure times > 7ks, to avoid observations 
that are too shallow to be useful for our purposes. We con- 
sider only sources with declination DEC(J2000) > — lOdeg 
to guarantee access from northern hemisphere telescopes 
and galactic latitude |bn| > 20deg to minimise contami- 
nation from Galactic stars. We use sources detected on the 
0.5-2 keV PN images (i.e. we do not consider the MOS1 and 
2 CCDs here) at off-axis angles < 14arcmin, to avoid spu- 
rious sources close to the edge of the field of view. 

T htt p : / / xmmssc- www . st ar . lc . ac . uk/ 

© 0000 RAS, MNRAS 000, 000-000 



Mining for normal galaxies in the 1XMM catalog 3 



Finally, in this paper we only consider fields with right 
ascensions, RA(J2000)> 4 hours, where our spectroscopic 
follow-up program of normal galaxy candidates is complete. 
For systems with RA(J2000) = — 4 hours optical spec- 
troscopy is in progress. A total of 51 XMM-Newton fields 
fulfil the criteria above. 

A drawback of the 1 XMM catalogue is that there is 
no areal coverage information providing an estimate of the 
surveyed area available at different point-source flux limits. 
To overpass this issue we determine the area curve indi- 
rectly by comparing the differential observed counts (i.e. in- 
cluding any incompleteness) in the 0.5-2 keV band with the 
expected logiV — log 5 from the literature. Here, we adopt 
the double-power law paramctrisation of Baldi et al. (2002) 
for the 0.5-2 keV differential source counts. For the observed 
counts we use all X-ray sources in the 1 XMM catalogue 
that overlap with the XMM-Newton fields considered in this 
paper and 0.5-2 keV likelihood probability > 7. The com- 
parison between our uncorrected counts and the Baldi et al. 
(2002) double-power law is shown in Figure 1. The result- 
ing area curve is plotted in Figure 2. This is increasing to 



uncorrected differential counts in Figure 1 are affected by (i) 
sources associated with the prime target of a given XMM- 
Newton pointing (Galactic stars or QSOs), (ii) small num- 
ber statistics and (iii) uncertainties arising from the extrap- 
olation of the Baldi et al. (2002) source counts to these 
fluxes. As a result the area curve above w 10 -13 ergs -1 cm~ 2 
is unreliable. However, this does not affect our results 
since all our normal galaxy candidates are fainter than 



10 



" 13 ergs 1 cm 



(see section 3) . We also confirm that the 
area curve in Figure 2 is robust using the Needles in the 
Haystack Survey (Georgantopoulos et al. 2005). For that 
sample we derive the area curve using the method outlined 
above and compare it with that estimated using the sensi- 
tivity maps. We find excellent agreement between the two 
methods suggesting that the area curve in Figure 2 is reli- 
able. 



2.2 Normal galaxy selection 

The SSC cross-correlates the 1 XMM X-ray source posi- 
tions with various astronomical catalogues including the 
U.S. Naval Observatory Catalog (USNO) version A2.0. This 
provides a homogeneous platform for the optical identifica- 
tion of 1 XMM sources. The USNO-A2.0 catalogue is based 
on photographic material and provides accurate positions 
(typical error of ~ 0.25 arcsec) and reasonable B-band pho- 
tometry with typical photometric errors of about 0.25 mag 
(Monet et al. 2003). This photometric uncertainty is ade- 
quate for this study that identifies 'normal' galaxy candi- 
dates on the basis of the low X-ray-to-optical flux ratio 
log fx I j opt < — 2. This is estimated from the 0.5-2 keV flux 
/x(0.5 - 2 keV) and the USNO-A2.0 B-band magnitude ac- 
cording to the relation 
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log /(0.5 - 2 keV) +0.4B + 5.29. 



(1) 



The equation above is derived from the X-ray-to-optical 
flux ratio definition of Stocke et al. (1991) that involved 
0.3-3.5 keV flux and V-band magnitude. These quantities 
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Figure 1. Differential source counts in the 0.5-2 keV band. The 
dots are the dN/ dS for all the 1 XMM catalogue sources with 
likelihood probability > 7 that lie on the XMM-Newton fields 
considered in this paper. The dashed line is the double power-law 
fit to the 0.5-2 keV counts of Baldi ct al. (2002). 
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Figure 2. Solid angle as a function of limiting flux in the 0.5- 
2 keV band. 



are converted to 0.5-2 keV flux and B-band magnitude as- 
suming a mean colour B — V = 0.8 and a power-law 
X-ray spectral energy distribution with index V — 1.8. 
The log fx / fopt — — 2 cutoff although minimises the AGN 
contamination may exclude from the sample X-ray ultra- 
luminous normal galaxies such as powerful starbursts (e.g. 
NGC3256; k 2 x 10 42 erg s _1 ; Moran, Lehnert & Helfand 
1999) or very massive ellipticals (Lx > 10 42 erg s _1 ; e.g. 
O'Sullivan, Forbes & Ponman 2001). Georgantopoulos et 
al. (2005) argue that in the case of star-forming galaxies 
the log fx I fopt = —2 limit is not likely to be a major 
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Figure 3. Cumulative normal galaxy counts in the 0.5-2 keV 
spectral band. Filled circles are the combined sample of nor- 
mal galaxy candidates from the Needles in the Haystack Survey 
(NHS; Georgantopoulos et al. 2005) and the present study. The 
open circles are the source counts of Hornschcmeier et al. (2003). 
The continuous line is the best fit to the log N — log S at bright 
fluxes from the combined NHS and 1XMM samples. The inlet 
plot shows the log TV — log S for the galaxies in the present study 
only, in comparison with the counts from the NHS 

source of incompleteness, at least for the low-z Universe. 
For early type galaxies we use the local E/S0 sample of Fab- 
biano et al. (1992) to estimate incompleteness because of the 
log fx /fopt = ~2 cutoff. We find that about 20 per cent of 
the galaxies in that sample have log Lx / Lb > —2, most of 
which are also luminous with Lx ii 10 42 ergs _1 . We note 
however, that some AGN contamination is expected within 
this X-ray bright subsample. 

The log fx /fopt < —2 regime is also populated with 
Galactic stars. The USNO A2.0 catalogue does not discrim- 
inate between optically extended and point-like sources. For 
this purpose we use the APM scans of the UKST plates* that 
provide star/galaxy separation. We further quantify the re- 
liability of the APM star/galaxy classification below, using 
follow-up spectroscopic observations or CCD quality data 
from the 3rd release of the Sloan Digital Sky Survey (SDSS 
DR3), where available. 

A total of 102 sources in our 1 XMM sample have 
log fx I fopt < —2. APM classifies 69 of them point-like and 
30 extended (e.g. galaxies) with the remaining 3 having am- 
biguous classification. The latter class is for objects that are 
assigned different types on the red and blue UKST survey 
plates. 

2.3 Optical spectroscopic observations 

A total of 21 out of 102 sources in our 1 XMM sample with 
log fx I fopt < —2 have been observed, 5 of them are unre- 

t http://www.ast.cam.ac.uk/~apmcat 



solved in APM with the remaining 16 classified galaxy-like. 
The 5 point-like sources are included in the target list for 
follow-up spectroscopic observations to explore the reliabil- 
ity of the APM star/galaxy separation. 

Low resolution optical spectroscopy was carried out 
with the Mexican 2-m class telescopes of the Observatorio 
Astrofisico Guillermo Haro (OAGH) in Cananea and Obser- 
vatorio Astronomico Nacional de San Pedro Martir (OAN 
SPM). Observations with the 2.1-m OAGH telescope were 
carried out with the Boiler & Chivens (B&Ch) spectrograph 
and Landessternwarte Faint Object Spectrograph and Cam- 
era (LFOSC) Zickgraf et al. (1997). 

The B&Ch spectrograph uses a Tektronix TK1024AB 
CCD mounted at the Cassegrain focus giving a pixel scale 
of approximately 0.45 arcsec. We use a 2.5 arcsec wide slit 
and a grating with 150 lines/mm providing a dispersion of 
3.5 A pixel" 1 and a wavelength resolution of ~ 15 A (~ 
4 pixels FWHM) over the range 3800-7100 A. The LFOSC 
is attached to the Cassegrain focus and is equipped with an 
EEV P8603 385x578 CCD giving a 10x6 arcmin 2 field of view 
and an image scale of 1 arcsec. We used a 3 arcsec wide slit 
and the G3 grism giving a wavelength resolution of ~ 18 A 
in the range 4200-9000 A. The observations were carried out 
during various observing runs between December 2003 and 
March 2005. 

The observations with the 2.1-m OAN-SPM telescope 
were carried out with the B&Ch spectrograph equipped with 
SITe3 (1024x1024 pix) CCD installed at the Cassegrain fo- 
cus giving a pixel scale of approximately 1 arcsec. We use a 
2.5 arcsec wide slit and a grating with 300 lines/mm, provid- 
ing a dispersion of 4.5A/pix and an effective instrumental 
spectral resolution of about 10 A (~ 2 pixels FWHM) in the 
wavelength range 3800-8000 A. The observations were car- 
ried out during two observing runs in November 2003 and 
February 2004. 

Typical exposure times for both the B&Ch and the 
LFOSC were 1 h per target. The observations were reduced 
following standard procedures using IRAF tasks resulting in 
wavelength and flux calibrated spectra. Redshifts were de- 
termined by visual inspection. Optical spectra of the galaxies 
from these observations are presented in Appendix 1. Our 
own follow-up spectroscopy is complemented with publicly 
available data from either the SDSS or the literature. 

Out of the 30 sources classified extended by APM one is 
the prime target of the XMM-Newton pointing (NGC 3184) 
and one has unresolved optical light profile in the SDSS 
CCD quality data and therefore is most likely associated 
with a Galactic star. Of the remaining 28 sources a total 
of 26 are assigned spectroscopic redshifts from either our 
own campaign (17) or the literature (9). Two sources in the 
sample have no redshifts. 

From the sources classified optically unresolved by the 
APM a total of 18 have either spectroscopic data from our 
own campaign or CCD quality star/galaxy separation from 
the SDSS. All these sources are indeed confirmed to be 
Galactic stars suggesting that the APM classification is re- 
liable. Finally, all 3 ambiguous sources are associated with 
Galactic stars on the basis of the optical spectra or SDSS 
photometric data. 
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3 THE SAMPLE 

In the analysis that follows we exclude all Galactic stars on 
the basis of either optical spectroscopy or the APM/SDSS 
star/galaxy separation. Our final normal galaxy sample 
comprises a total of 28 sources, 26 spectroscopically iden- 
tified galaxies and 2 sources without spectra classified ex- 
tended by APM. This is presented in Table 1. For complete- 
ness we also show the 4 sources with extended (1) or ambigu- 
ous (3) APM classification that turned out to be Galactic 
stars. 

A number of sources in our sample are classified AGNs 
on the basis of their optical spectroscopic properties. These 
include source #A063 which is a BAL-QSO at z = 0.149 
(Gallagher et al. 1999) and source #AI40 (NGC 4156) sug- 
gested to harbor AGN activity by Elvis et al. (1981). Two 
of the systems for which we obtained optical spectroscopic 
observations, #A001 and A035 (see Fig. Af), show evi- 
dence for broad Ha emission-line with FWHM Sv « 1600 
and 2000 km s -1 respectively. In both cases the measured 
broad width may be partly due to the low resolution spec- 
trum that does not allow separation of the Ha from the 
[Nil] 6583 A. Nevertheless, we conservatively classify these 
sources AGN. Source #A019 has narrow optical emission 
lines but the line ratios [log ([S II] 6716 + 31/Hq) « -0.4, 
log([OIII]5007/H/3) « +0.5; sec Fig. Al] place it in the 
AGN region of the diagnostic diagrams of Kewley et al. 
(2001) first introduced by Baldwin, Phillips & Terlevich 
(1981) and Veilleux & Osterbrock (1987). Here we adopt 
the theoretical lower bound for starburst galaxy emission- 
line ratios from Kewley et al. (2001; see their Figure 16) to 
discriminate between H II and AGN dominated systems. All 
these sources are excluded from the analysis. Additionally 
two of our normal galaxies (sources #A106, A149; see Fig. 
Al) are Coma cluster members identified in fields targeting 
this cluster and are excluded from statistical studies (e.g. 
log N — log S, luminosity function). 

Table 1 presents our sample including AGNs. We list: 

1. Identification number and the NED name of that 
source if available. 

2, 3. Right ascension and declination of the optical 
source in J2000. 

4. B-band magnitude from the USNO-2 version A2.0. 

5. if-band magnitude from the 2MASS All Sky Data 
Release (Skrutskie et al. 1997). Either the extended or the 
point source catalogue was used. 

6. Offset in arcseconds between the X-ray and optical 
source position. 

7. 0.5-2 keV X-ray flux corrected for Galactic absorption 
in units of 10 -14 ergs -1 cm~ 2 . 

8. Hardness ratio using the 0.5-2 keV and the 2-4.5 keV 
spectral bands corrected for vignetting. 

9. X-ray-to-optical flux ratio defined in equation 1. 

10. 1.4 GHz radio flux density in mjy from either the 
FIRST (Becker et al. 1995; White et al. 1997) or the NVSS 
(Condon et al. 1998) surveys. 

11. Redshift of the source. In the appendix we present 
the optical spectra of the normal galaxy candidates in the 
1XMM sample obtained at the OAGH and OAN-SPM Mex- 
ican telescopes as part of this project. 

12. 0.5-2 keV X-ray luminosity in units of ergs -1 . For 



the k-correction a power-law spectral energy distribution 
was adopted with photon index F = 1.8. 

13. Ha luminosity in units of ergs -1 . This is measured 
from the optical spectra after correcting for intrinsic dust 
obscuration using the Balmer decrement Ha/H/3. For more 
details see section 6. 

14. 1.4 GHz radio luminosity. For the k-correction a 
power-law spectral energy distribution was adopted with 
spectral index a = 0.8. 

15. Classification on the basis of the optical spectro- 
scopic observations: HII for starforming galaxies, ABS for 
absorption line spectra and AGN for systems showing evi- 
dence for central black hole accretion. There is no spectrum 
available for source #142 with the redshift estimate, show- 
ing both emission and absorption line features, coming from 
Arp (1977). In the analysis that follows this is assumed to 
be H II type system. 

16. APM star/galaxy separation. 



4 THE log N — log S 

The inlet plot in Figure 3 presents the log N — log S in the 
0.5-2 keV band for our normal galaxy candidates derived 
using the area curve from section 2.1. The 5 AGNs and 
the two Coma cluster members in Table 1 are not used 
in this plot. At faint fluxes (<, 10 -14 erg s -1 cm -2 ) there 
is fair agreement with the results from the Needles in the 
Haystack Survey (NHS; Georgakakis et al. 2004; Georgan- 
topoulos et al. 2005) that used normal galaxies selected in 
the 0.5-8 keV band. The NHS cumulative counts are shifted 
to the 0.5-2 keV energy range assuming F = 1.8. At the 
bright end our counts are elevated compared to this survey, 
although still consistent within the la uncertainties. This 
may be because of field-to-field variations affecting primar- 
ily the bright flux counts derived from our relatively small 
area survey (« 6deg 2 ). Despite these uncertainties at the 
bright end, the agreement of our log N — log S with previous 
studies suggests that the area curve derived in section 2.1 is 
robust. 

We improve the statistical reliability of our results by 
merging our normal galaxy sample with that from the NHS 
presented by Georgantopoulos et al. (2005) . As already men- 
tioned this sample is compiled in the 0.5-8 keV spectral band 
using public XMM-Newton fields overlapping with the SDSS 
DR-2. A total of 28 systems over « lldeg 2 are identified. 
We shift the NHS 0.5-8 keV fluxes and area curve into the 
0.5-2 keV band adopting T = 1.8. A total of 13 overlapping 
fields in the two surveys are taken into account only once 
in the area curve calculation. The combined sample covers 
about ~ 15 deg 2 and comprises a total of 46 unique entries. 
The logiV — logS is presented in Figure 3. The differential 
'normal' galaxy counts in the 0.5-2 keV of the two combined 
surveys are fit by a power law yielding a slope of — 2.46±0.13. 
This correspond to a slope of — 1.46 ± 0.13 for the cumula- 
tive number counts, in agreement with the results of Tajer 
et al. (2005) at bright fluxes and Hornschemeier et al. (2003; 
— 1.46t^3o) at the faint end, the latter sample based on the 
same selection criteria used here (log/x// op < < —2). 
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The columns are: 1: identification number and common galaxy name if available; 2: right ascension of the optical source; 3: declination of the optical source; 4: optical S-band magnitude; 
5: 2MASS K-band magnitude; 6: offset in arcsec between the optical and X-ray source position; 7: X-ray flux in the 0.5-2 keV spectral band in units of 10 -14 erg s -1 cm -2 ; 8: Hardness ratio 
estimated from the 0.5-2.0 and 2.0-4.5 keV bands ; 9: X-ray-to-optical flux ratio; 10: radio flux density; 11: spectroscopic redshift; 12: 0.5-2 keV X-ray luminosity in units of ergs -1 ; 13: Ha 
luminosity in units of ergs -1 ; 14: 1.4GHz radio power in units ofW/Hz; 15: classification based on optical spectroscopic properties; 16: APM star/galaxy separation on the basis of the optical 
light profile. 

a There is no spectrum available for source #142 with the redshift estimate, showing both emission and absorption line features, coming from Arp (1977). In the analysis that follows this is 
assumed to be HII type system. 



Table 1: The candidate 'normal' galaxy sample. 
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5 THE LUMINOSITY FUNCTION 

We further explore the statistical properties of the normal 
galaxy sample by deriving the X-ray luminosity function 
(LF) using methods that are fully described by Georgan- 
topoulos et al. (2005). We estimate both the binned X-ray 
LF using the technique of Page & Carrera (2000) and the 
parametric Maximum Likelihood fit (Tammann, Sandage & 
Yahil 1979) adopting a Schechter (1976) form for the lumi- 
nosity function. 

We first estimate the 0.5-2 keV LF for the IXMM data 
alone. Figure 4 plots our results in comparison with those 
from the combined NHS and CDF data from Georgantopou- 
los et al. (2005) shifted to the 0.5-2 keV band assuming 
r = 1.8. There is good agreement within the la errors fur- 
ther suggesting that the area curve derived in section 2.1 is 
adequate for statistical studies. We note that the two sources 
without spectroscopic redshift estimate are not used in this 
calculation. 

We improve the statistical reliability of our LF esti- 
mates by combining our sample with both the NHS and 23 
z <J 0.2 galaxies selected in the 0.5-2 keV spectral band from 
the CDF-N and CDF-S. The CDF-N data are obtained from 
Hornschemeier et al. (2003) by selecting a total of 15 sources 
with 0.5-2 keV band detection. All of these systems have 
spectroscopic redshifts available. In the case of the CDF-S 
we select a total of 8 sources detected in the 0.5-2 keV spec- 
tral band with log fx / f opt < —2 from the catalogue pre- 
sented by Rosati et al. (2001). Spectroscopic (total of 5) or 
photometric (total of 3) redshifts are available from Szokoly 
et al. (2004) and Zheng et al. (2004) respectively. This com- 
bined sample is also split into absorption and emission line 
systems to explore the LF for different galaxy types. Our 
results are plotted in Figure 5 with the maximum likelihood 
best-fit parameters presented in Table 2, where L* is the 
break of Schechter function, a is the faint-end slope and 
<f>-t, the normalisation. In the same table we give the X-ray 
emissivity (luminosity per Mpc ) estimated by the relation 
jx = / &(L) L dL, where $(L) is the luminosity function. 
The la errors on and a are estimated from the regions 
around the best fit where the likelihood function changes by 
SL — 0.5 (e.g. Press et al. 1992). The uncertainties in 0* and 
jx are approximated by performing 200 bootstrap resamples 
of the data and then estimating the 68th percentiles around 
the median. For a Gaussian distribution these correspond to 
the 68 per cent confidence level. 

Figure 5 also compares the X-ray luminosity function 
of the combined 1XMM+NHS+CDF sample at a mean red- 
shift z = 0.087 with the higher- 2 LF estimates of Norman 
et al. (2004). These authors used Bayesian statistical anal- 
ysis to select normal galaxy candidates in the CDF-North 
and South and derived the first ever X-ray luminosity func- 
tion for these systems in two redshift bins with medians 
z — 0.26 and z = 0.66 respectively. Their z < 0.5 LF in 
Figure 5 is in fair agreement with ours at the faint end but 
diverges at brighter luminosities. This may suggest (i) con- 
tamination of the Norman et al. (2004) sample by AGNs 
at bright luminosities, (ii) bias in our sample against X-ray 
ultra- luminous systems (e.g. Lx ii 10 42 ergs _1 ) because of 
the log(/x// op t) < —2 selection or (iii) evolution of the 
normal galaxy luminosity function. In the latter case, the 
median redshift of the z < 0.5 subsample of Norman et al. 




39 40 41 42 

log L x (0.5-2 keV) (erg s" 1 ) 

Figure 4. The local 0.5-2 keV luminosity function for the present 
sample of IXMM galaxies. Also shown is the maximum likelihood 
fit method for z < 0.2 normal galaxies from the combined NHS 
and CDF samples (see Georgantopoulos et al. 2005) shifted from 
the 0.5-8 to the 0.5-2 keV band assuming T = 1.8. 



is Zmedian = 0.26 higher than our mean of z — 0.087. For lu- 
minosity evolution of the form (1 + z) 2 ' 7 derived by Norman 
et al. (2004), a source at z = 0.26 is expected to become 
about 1.5 times more luminous relative to z = 0.087. Such 
a brightening can indeed, partially account for the observed 
differences. 

In Table 2 we also estimate a 0.5-2 keV X-ray emissivity 
of (0.24 ± 0.02) x 10 38 and (0.22 ± 0.03) x 10 38 ergs" 1 Mpc" 3 
for emission and absorption line systems respectively. These 
values are lower than previous estimates. For example Geor- 
gakakis et al. (2003b) used stacking analysis to estimate 
the mean X-ray properties of 2dF galaxies at z ~ 0.1. 
These authors estimate j x = (0.4 ± 0.3) x 10 38 and (1.2 ± 
0.4) x 10 38 ergs -1 Mpc -3 for late and early type galaxies 
respectively. Georgantopoulos, Basilakos & Plionis (1999) 
convolved the local optical luminosity function of the Ho, 
Filippenko & Sargent (1997) sample with the correspond- 
ing Lx — Lb relation based on Einstein data. Their jx es- 
timates scaled to the 0.5-2 keV band are (0.50 ± 0.06) x 
10 38 ergsec -1 Mpc -3 for HII galaxies and (0.53 ± 0.06) x 
10 38 ergsec _1 Mpc~ 3 for passive galaxies. The origin of the 
discrepancy between the emissivity estimates presented here 
and those found in the studies above is not clear. In the case 
of stacking analysis for example, luminous systems may bias 
the mean signal toward high values, resulting in an overes- 
timation of the jx ■ Although the method presented here for 
estimating jx is direct and robust, we note that a possible 
bias in our sample against high Lx systems may result in a 
systematic reduction of the jx ■ 
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Sample Number of mean log Z,* a <f>+ j x 

sources redshift (ergs" 1 ) (xln(10) X 10~ 4 Mpc~ 3 dcx" 1 ) (xl0 38 crgs^ 1 Mpc" 3 ) 

1XMM/NHS/CDF 67* 0.087 41.02to.la - 1 - 76 -°oao 3 - 4 ±i'.8 

Emission 33 0.091 -1.61^° e -°t 7 3 .l 

Absorption 34 0.083 41.25±°;f| -1. 791^4 °- 9 -o.i 

*The 2 sources in the 1XMM sample without redshift determination have been excluded from the analysis. 



0.47 
0.24 
0.22 



+0.04 
0.04 
+0.02 
0.02 
+0.03 
0.03 



Table 2. The luminosity function best-fit parameters 




39 40 41 42 43 



log L x (0.5-2 keV) (erg s" 1 ) 

Figure 5. The 0.5-2 kcV luminosity function for the combined 
sample of 1XMM, NHS and CDF z < 0.2 normal galaxies. The 
derived maximum-likelihood fits to the observations arc shown 
with the continuous (full sample), dotted (emission-line objects) 
and dashed (absorption-line systems) lines. For clarity the results 
from the non-parametric method are shown with filled circles only 
for the full sample. The luminosity function derived by Norman 
ct al. (2004) in two redshift bins is plotted for comparison. 

6 X-RAY LUMINOSITY AND HOST GALAXY 
PROPERTIES 

Using the 1XMM+NHS sample we explore the correlation 
between 0.5-2 keV X-ray luminosity and indicators of star- 
formation activity. The Ha luminosity (Lho) measured from 
the optical spectra and the 1.4 GHz radio luminosity (£1.4) 
obtained from either the FIRST or the NVSS surveys are 
used as proxies to the galaxy SFR. Only X-ray selected nor- 
mal galaxy candidates with emission line optical spectra are 
considered here. 

Figure 6 plots Lx against L1.4 for our sources. Because 
of the large number of upper limits we do not attempt to fit 
an Lx — L1.4 relation to the data. Also shown in this figure 
for comparison are the spiral galaxies from Shapley et al. 
(2001). These authors compiled an optically selected sample 
of local spirals with X-ray data (detections or upper lim- 
its) from the Einstein observatory and 4.85 GHz radio flux 
densities from the literature. Powerful AGNs in this sam- 
ple are not plotted in Figure 6. The Einstein X-ray fluxes in 
the 0.2-4 keV band are transformed to the 0.5-2 keV band as- 
suming a power-law model with a spectral index T = 1.8. At 



radio wavelengths a power-law spectral energy distribution 
of the form f„ ~ i/~ 0,8 , appropriate for star- forming galax- 
ies, is adopted to convert the 4.85 GHz radio flux density of 
Shapley et al. (2001) to 1.4 GHz. Also shown in this plot is 
the Lx — L\a relation of Ranalli et al. (2003) of the form 
Lx oc Li l 8 . These authors used ASCA and BeppoSAX data 
of optically selected star-forming galaxies classified on the 
basis of high quality nuclear spectra from Ho et al. (1997). 
Many of our sources are not detected at 1.4 GHz to the 1 mjy 
flux density limit of the FIRST survey. The upper limits 
and detections, although within the scatter of the Shapley 
ct al. (2001) data, appear slightly offset from the Ranalli 
ct al. (2003) best-fit relation, particularly for star-formation 
rates < 10 MQ yr _1 . This may suggest a contribution from 
low-mass X-ray binaries (e.g. old stellar population) to the 
observed X-ray emission of more quiescent systems. A larger 
number of radio detections is required to further explore this 
issue. 

In Figure 7 we plot Lx against Ha luminosity for 
emission-line galaxies in the combined 1XMM+NHS sample 
only. The latter is measured from the optical spectra after 
correcting for intrinsic dust obscuration. We quantify the 
visual extinction, Ay, by comparing the Ha/H/3 decrement 
with the theoretical (case B recombination) value of 2.86 
(Brocklehurst 1971) and a standard reddening curve (Sav- 
age & Mathis 1979). The technique is hampered by the poor 
S/N ratio of some of the spectra and by stellar H/3 absorp- 
tion that reduces the measured flux of this line. We account 
for this effect by applying a correction of —2 A to H/3, sim- 
ilar to the mean stellar absorption in this line determined 
for star-forming galaxies (Tresse et al. 1996; Georgakakis et 
al. 1999). A small number of systems in our sample have 
Ha/H/3 ratios below the theoretical value of 2.86 for case B 
recombination, suggesting either little dust reddening, stel- 
lar H/3 absorption lower than —2 A or poor signal-to-noise 
ratio optical spectra. For these sources we still apply a cor- 
rection to Ha for dust obscuration adopting the mean optical 
extinction of Av = 1 mag for spiral galaxies determined by 
Kennicutt (1992). 

The best-fit relation of the form Lx oc L^ a with 
/3 = 0.69 ± 0.06 is plotted in Figure 7. This is estimated 
using the bisector least square fitting described in Isobe et 
al. (1990). Also shown is the Lx — Ln a relation for star- 
forming galaxies determined by Zezas (2000) using a total 
of 43 systems with PSPC ROSAT data, classified as H II 
on the basis of high quality nuclear optical spectra from Ho 
et al. (1997). The exponent of the Zezas (2000) relation is 
/3 = 0.62 ±0.11 in good agreement with our determination. 
There is however a difference in the normalisation. This is 
most likely because Zezas (2000) uses Ha luminosities inte- 
grated over the whole galaxy whereas we are measuring this 
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Figure 6. 0.5-2 keV X-ray luminosity against radio power. The 
L1.4 is converted into SFR on the top horizontal axis using the 
relation of Bell ct al. (2003). The large open squares are the 
emission-line galaxies from the 1XMM+NHS sample. A cross on 
top of a symbol is for AGNs. Radio upper limits are shown with 
an arrow pointing to the left. The small open circles are local 
spiral galaxies from Shaplcy ct al. (2001). The dashed line is the 
best fit Lx — Li.4 relation of the form Lj oc L^' 88 derived by 
Ranalli ct al. (2003). The dotted lines represent the 1 sigma rms 
envelope around the best fit. 

line through a slit (or a fibre in the case of SDSS data) with 
typical angular size of 2 arcsec. This width corresponds to a 
physical scale of about 2 kpc diameter at the mean redshift 
of our sample. Extrapolating fluxes measured through a slit 
to integrated values over the whole galaxy is not straightfor- 
ward and depends on the size of star-forming regions, their 
distribution in the galaxy, the position of the slit and the 
seeing at the time of the observations. 

We further explore the association between 0.5-2 keV X- 
ray luminosity and galaxy mass for absorption-line galaxies 
in the combined 1XMM+NHS sample. The 7v-band lumi- 
nosity is used as proxy to the stellar mass of the system. 
Figure 8 plots Lx against Lk- Also shown in this figure 
is the relation between the total X-ray luminosity of point 
sources, Lxp, and the host galaxy Lk determined by Col- 
bert et al. (2004) using Chandra observations of nearby el- 
lipticals. This line shows the mean expected contribution of 
low-mass X-ray binaries to the observed emission at a given 
Tf-band luminosity or galaxy mass and does not include the 
X-ray emitting hot gas component. 

For absorption- line systems in our sample (excluding 
AGNs) we adopt a relation of the form Lx oc L K to esti- 
mate a bisector exponent 7 = 1.5 ± 0.1. This is flatter than 
the slope of « 1.8 found by Fabbiano et al. (1992) for the 
Lx — Lb relation of local E/SOs and Hornschemeier et al. 
(2005) for the Lx /stellar- mass correlation of absorption line 
galaxies in the SDSS. This apparent discrepancy may be due 
to the log fx I j opt < —2 cutoff that effectively excludes from 
the sample very luminous X-ray selected normal galaxies. 
For example about 20 per cent of the Fabbiano et al. (1992) 
sample have log Lx/Lb > —2. The vast majority of these 



Figure 7. 0.5-2 keV X-ray luminosity against Ho luminos- 
ity corrected for intrinsic dust obscuration for the combined 
1XMM+NHS emission-line galaxies. The Ho is converted into 
SFR on the top horizontal axis using the relation of Kcnnicutt 
(1998). A cross on top of a symbol is for AGNs. The continu- 
ous line shows the bisector least squares fit relation of the form 
L x oc L f ^ a with (3 = 0.69 ±0.06. The dotted line is the L x - L Ha 
relation for star- forming galaxies of Zczas (2000). Although the 
slopes of the two lines are similar they differ in the normalisation. 
This is most likely because Zczas (2000) used Ha data integrated 
over the whole galaxy rather than slit spectroscopy used here. 

systems are also luminous with Lx ii lO^ergs^ 1 . Exclud- 
ing E/S0 galaxies with log Lx /Lb > —2 from that sample 
we find a flatter slope of « 1.5 for the Lx — Lb relation 
in agreement with the value estimated here. Alternative the 
flatter Lx — Lk relation may indicate that the if-band pro- 
vides a better proxy to galaxy mass. For example, Shapley 
et al. (2001) used ff-band near-infrared data and find a flat- 
ter slope for the Lx — Lh relation of spirals compared to the 
Lx ~Lb relation for the same systems. Although comparing 
the X-ray properties of ellipticals and spirals is not appro- 
priate it suggests that the flatter slope we are estimating 
may be partly due to the use of near-infrared rather than 
optical luminosities. 



7 SUMMARY AND CONCLUSIONS 

In this paper we demonstrate the power of the First XMM- 
Newton Serendipitous Source Catalog for studies of X-ray 
selected normal galaxies. Our sample is compiled from X- 
ray sources detected on XMM-Newton pointings that have 
(i) EPIC-PN detector as prime instrument operated in full- 
frame mode, (ii) exposure time > 7ks, (iii) declinations 
DEC(J2000) > -lOdeg, (iv) galactic latitude |b n | > 20deg 
and (v) right ascension RA(J2000) > 4 hours. A total of 
51 fields fulfil the above criteria, covering a total area of 
» 6deg 2 to the 0.5-2 keV band limit /x(0.5 - 2keV) « 
lO^ergs^cm- 2 . 

The USNO A2.0 catalogue provides a homogeneous 
platform that allows identification of X-ray selected normal 
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Figure 8. 0.5-2 kcV X-ray luminosity against A"-band luminos- 
ity. The large squares are normal galaxies in the NHS and 1XMM 
samples with absorption line optical spectra. A cross on top of a 
symbol is for sources in our sample that show evidence for AGN 
activity. The continuous line is the Lx — Lk relation of Colbert 
ct al. (2004) which, for a given if-band luminosity, corresponds 
to the mean expected contribution of the point source X-ray bi- 
nary population to Lx ■ The dashed line is the best fit relation to 
absorption line galaxies in our combined NHS+1XMM sample. 



galaxy candidates on the basis of the low X-ray-to-optical 
flux ratios of these sources, log fx / fopt < —2. Reliable 
star/galaxy separation is available from the APM, allowing 
us to exclude Galactic stars that also have log fx/fo P t < — 2. 
Optical spectroscopy is then used to identify systems that 
show evidence for AGN activity resulting in a sample of 
23 normal galaxy candidates: 9 with narrow emission lines, 
12 with absorption lines only and 2 without optical spec- 
troscopic information. Future releases of the XMM-Newton 
Serendipitous Source Catalog will provide much wider areal 
coverage resulting in significantly larger \ow-z normal galaxy 
samples. 

At present we increase our sample size by combin- 
ing it with X-ray selected normal galaxy candidates from 
the Needles in the Haystack Survey (Georgantopoulos et 
al. 2005). This provides a total of 46 z <, 0.2 X-ray de- 
tected normal galaxies, the largest low-z sample yet avail- 
able. Such a large number of sources provides a unique op- 
portunity to constrain the normal galaxy log N — log S at 
bright fluxes (10~ 15 — 10~ 13 erg s _1 cm -2 ). We estimate a 
slope of —1.46 ± 0.13 consistent with the euclidean predic- 
tion and in agreement with previous determinations (Horn- 
schemeier et al. 2003; Tajer et al. 2005). 

Our enlarged sample is further combined with 23 local 
(z ^ 0.2) galaxies from the Chandra Deep Field North and 
South surveys to construct the local X-ray luminosity func- 
tion of normal galaxies. This is fit with a Schechter function 
with a break at logL* = 41.02^q ergs -1 and a slope of 
a = —1.76 ± 0.10. The large sample size allows us to also 
estimate the luminosity function of emission and absorption 
line systems separately and to provide the most accurate es- 



timates of the 0.5-2 keV X-ray emissivity of these systems in 
the nearby Universe. We find j x = (0.24 ± 0.02) x 10 38 and 
(0.22 ± 0.03) x 10 38 erg s -1 Mpc~ 3 for emission and absorp- 
tion line systems respectively. We note that these are lower 
than previous (less direct) estimates based on stacking anal- 
ysis (e.g. Georgakakis et al. 2003b) or optically selected local 
galaxies (Georgantopoulos et al. 1999). 

Finally, for the combined sample of 46 systems we ex- 
plore the association between X-ray luminosity and host 
galaxy properties, such as SFR and stellar mass. We use 
the 1.4 GHz radio and Ha luminosities as star- formation in- 
dicators. We find that the Lx — Li.4 relation for our sources 
is consistent with that determined by Ranalli et al. (2003) 
for the local Universe. For the log Lx — log Lua correlation 
we estimate a slope of ~ 0.7 in fair agreement with previous 
results (Zezas 2000). The above relations suggest that the 
X-ray luminosity in the emission-line galaxies in our sample 
is directly related to the current star-formation activity as 
measured by the 1.4 GHz radio and Ha luminosities. For the 
early type galaxies in the sample showing absorption optical 
lines only we use the Jf-band as proxy to their stellar mass 
and find a linear correlation between log Lx and log Lk with 
a slope of ~ 1.5. This is flatter than the slope of ~ 1.8 for the 
Lx—Lb relation for local ellipticals of Fabbiano et al. (1992) 
or the Lx / stellar-mass correlation for absorption line galax- 
ies in the SDSS determined by Hornschemeier et al. (2005). 
This may be due to a possible bias in our sample against 
very luminous galaxies, Lx > 10 42 ergs _1 cm~ 2 , likely in- 
troduced by the low X-ray-to-optical flux ratio selection. 
Alternatively, this may indicate that the 7\-band provides a 
better representation of the galaxy mass. 
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APPENDIX A: 

In this appendix we present the optical spectra of the nor- 
mal galaxy candidates in the 1XMM sample obtained at the 
OAGH and OAN-SPM Mexican telescopes as part of this 
project. 
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Figure Al. Optical spectra of galaxy candidates in the 1XMM sample obtained at the OAGH and OAN-SPM Mexican telescopes. 
Sources A001, A035 and A019 are classified AGN. 
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Figure Al - continued 
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